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COSMIC RAY ISOTOPES* 

E. C. Stone 

i California Institute of Technology 

Pasadena, California (USA) 91109 


Since this is the first rapporteur talk on cosmic ray iso- 
topes in thirteen Cosmic Ray Conferences, I will begin by 
briefly summarizing our interest in isotopic composition 
as follows : 

a) Isotopes, not elements, are the result of various 
nucleosynthesis processes. Thus, we must know the 
isotopic composition of cosmic rays in order to 
adequately exploit the relationship of cosmic rays 
and stellar processes. 

. ‘ \ 

b) Isotopes, not elements, are the result of second- 
ary fragmentation processes, either in cosmic ray 
source regions or in the interstellar medium. 

Thus, we must know the isotopic composition in 
order to completely unravel source composition 
from secondary composition. 

c) Isotopes, not elements, are radioactive. Thus, 
we must know the isotopic composition in order to 
take full advantage of various radioactive nuclei 
as cosmic ray clocks which can be used to determine 
the age of cosmic rays. 

With this short summary, I now turn to a brief review of 
the results presented at this conference. Since I can’t 
possibly discuss each new measurement, I have chosen the 
following topics: I) H and He isotopes, II) Be isotopes, 

III) A1 isotopes, IV) Fe isotopes, V) cross sections and 
model calculations, and VI) instrumental considerations. 


I. 

H and He 

Isotopes. The 

relevant results are reported in the 

follow- 

ing 

papers : 





Session 

Paper 

Authors 



0G-3A 

400 

Teagarden, von Rosenvinge, McDonald 

- GSFC-1 



159 

Hurford, Mewaldt, Stone, Vogt 

- CIT 



375 

Apparao 

- Tata 


^Rapporteur paper, presented at the 13th International Cosmic Ray Conference, 
Denver, August 26, 1973. 
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Study of the H and He isotopes can provide information on the interstellar 
spectra of H and 4 He, because D and 3 He are thought to be secondaries pro- 
duced almost entirely by galactic cosmic ray and ^He. In particular, 
the ratio of D/^He is especially sensitive to the primary spectra because 
low energy D is copiously produced by high energy protons through spall- 
ation of interstellar ^He and by the reaction p+p -*■ D-Hr. Thus, the ratio 
D/'+He at’ low kinetic energies is really a ratio of a high energy flux (the 
which produced the D) to a low energy flux (the 4 He) and is therefore 
sensitive to the interstellar spectrum. 


At the Hobart Conference, J. P. 
Meyer (1971) compared extensive 
calculations of D/^He with data 
and concluded that the inter- 
stellar spectra were close to 
power laws in total energy. The 
data available at that time 
covered the period from - 1967 
to 1970. 

At this conference strikingly 
different results have been 
obtained for 1972 by the Caltech 
and Goddard groups. These new 
results are shown in Figure lj 
along with the earlier results 
summarized by Simpson (1971). 

Several points are contained in 
Figure 1. First, the ratio 
D/^He was much smaller in 1972 
than in previous years. This 
reduced ratio is partly due to 
a factor of three reduction in 
the D flux compared to 1967 
and partly due to an enhanced 
4 He flux which is larger than 
that observed at the last solar 
minimum. 



Fig, 7 . The deuterium to ^He ratio. 
The new 1972 results are 
compared with previously 
summarized results and 
calculations. 


The second key point is that we expect adiabatic deceleration to reduce the 
energy dependence of the ratio below ~ 200 MeV/nucleon (see, for example, 
Biswas and Ramadurai, paper 406). Adiabatic deceleration occurs because 
as the galactic particles penetrate into 1 AU they are diffusing through 
the expanding interplanetary medium, which produces a "cooling" effect 
.similar to the cooling of a gas in an expanding box. Thus, the 10 to 30 
MeV/nucleon particles at 1 AU would all have had higher energies (> 100 
MeV/nucleon) outside the modulation region and the ratio D/'+He below 100 
MeV/nucleon should be very similar to that at ~ 200 to 500 MeV/nucleon. 




-3- 


Figure 1 includes the interstellar ratios calculated by Meyer (1971) for 
power laws of the form (T+T ) -2,6 for the primary 1 H and 4 He spectra, 

where T is kinetic energy and T q =0 and 500 MeV. The predicted ratio is 

~ 0.1 in the 200 MeV/nucleon region for both assumed power laws. Thus, 
adiabatic deceleration should have resulted in a similar ratio of ~ 0.1 
at lower energies, which is in serious disagreement with the 1972 observa- 
tions. 

One possible explanation is that the propagation of <50 MeV/nucleon par- 
ticles is not dominated by adiabatic deceleration, but that some fraction 
of both low energy D and 4 He are gaining access without large energy losses 
at times near solar minimum. Should this be the case, then the energy 
dependence of the ratio at 1 AU might directly reflect the unmodulated 
energy dependence. As shown in Figure 1, there is a striking similarity 
between the calculated interstellar ratio for the kinetic energy power law 

case (T =0) and the observations at 1 AU. 
o 

Another possibility is that either the interstellar D or 4 He spectra is 
much different than assumed in the model calculations. Perhaps a signif- 
icant enhancement in the low-energy interstellar 4 He could, even with 
adiabatic deceleration, produce a noticeable increase in 4 He at 1 AU with- 
out a corresponding increase in the production of D. 

The proper interpretation of these very recent results may be quite differ- 
ent once the appropriate calculations have been made. There is, however, 
a strong possibility that our understanding of solar modulation may be 
significantly modified and that we may obtain important information on the 
interstellar spectra of low energy H and 4 He. 

To close off this subject, I want to mention that Apparao reported an upper 
limit of D/He< 0.6 above 16.8 GV. In the revised version of the paper, 
the author points out that the proton background seriously limited the 
results. An improvement of a factor of 5 to 10 in sensitivity would be 
of considerable interest, since the D/He ratio at higher energies depends 
on the proton pathlength in the interstellar medium. ' 

IX. Be Isotopes. The relevant results are reported in the following 
papers : 

Session Paper Authors 

OG-3A 156 Garcia-Munoz , Mason, Simpson (Chicago) 

329 Enge, Fukui, Beaujean (Kiel) 
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The primary interest here, of course* is in the 10 Be isotope, which has 
a radioactive halflife of ~ 1.5 x 10 6 y and thus is an appropriate cosmic 
ray clock. However, 10 Be is difficult to separate from *Be, and fluxes 
of only 7 Be and 9 Be were resolved by the Chicago group. The Kiel group 
gave a status report of their investigation which is not yet complete. 


The Chicago results for Be are shown 
in Figure 2, which illustrates the 
mass resolution obtained. Note that 
the separation of 7 Be and 9 Be is 
similar to that of ?He and ^He. 

There is however, no statistically 
significant indication of a 10 Be 
peak, and the authors emphasize only 
that they have clearly separated 7 Be 
and 3 Be. The authors do report 3 
events in the 10 Be region in excess 
of that expected from 9 Be contam- 
ination and background. 

The Chicago results are 

7 Be 

~~ = 0.50+0.07 


^Be 

Be 


0.4110.10 


10 Be 

= 0 .09±0 . 10 


<°) 0« T B.* Be 0 



c dE/dx Channel number 



dE/dx Channel number 


which are consistent with galactic 
cosmic ray propagation calculations. 

Because of the large uncertainty of 
the 1( ^Be flux, no statement can be 
made about cosmic ray lifetime, 

Although not reported at this conference, I do feel that it is appropriate to 

mention recent results by Webber, Lezniak, Kish, Damle (1973) obtained with 

a balloon borne instrument. Their results show separation of g Be and 10 Be 

and they report 7 1( ^Be nuclei out of 100 Be nuclei observed at 2.9 g/cm 2 

atmospheric depth. When adjusted for equal energy intervals and for solar 

modulation effects, they calculate that they should have observed 18±3 10 Be 

nuclei if none had decayed. Thus, they conclude that the mean lifetime 

T of cosmic rays is 
cr 

T - 3.4±?*t x 10 6 y 
cr 1.3 ■ 


Fig, 2. The observed mass 

distributions for Be 
and He reported by the 
Chiaago group . 
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This is probably the first indication that we will indeed be able to use 
i0 Be as a clock with further improvements in statistical accuracy, frag- 
mentation parameters, and instrumentation. It is especially important 
to accumulate at least ~ 100 10 Be nuclei with a given instrument so that 
any residual instrumental pecularities or background can be quantified 
and eliminated. 

III. A1 Isotopes. The relevant papers are the following: 

Session Paper Authors 

0G-3A 305 Beaujean, Enge (Kiel) 

043 Webber, Lezniak, Kish (UNH) 

Another possible cosmic ray clock is 26 A1, which has a halflife of ~ 7xl0 5 y. 
Since secondary 2 ^A1 is the result of the interaction of heavier primaries 
(e.g.,Fe, Ca, Si), while 10 Be is produced mainly by lighter primaries (C,0), 
comparison of these two clocks might indicate whether or not the residence 
time for the two groups of primaries is the same. 

The Kiel group described their work with plastic track detectors which may 
ultimately yield a mass resolution of a ~ 0.3 amu. Presently the resolution 
is ~ 1 amu and no comment is possible on the 2 ®A1 abundance. 

The UNH group reported results on the isotopic abundances of a broad range 
of elements. For A1 they report 

2S A1+ 26 A1 = 2 
27 A1 " 9 

where the numbers in the ratio correspond to the observed number of events. 
They calculate that secondary production would yield at ratio of 26 A1/ 27 A1 
^ 0,35* The observed ratio would be smaller due to either the decay of 
26 A1 or the presence of 27 A1 in the source. Thus the observed ratio is 
suggestive of the presence of 26 A1, but the uncertainties in the data are 
presently too large to support any conclusion about the decay of 26 A1* 

The observed mass distribution for A1 is included in Figure 3 in the next, 
section* 

IV. Fe Isotopes* The relevant papers are: 

Session Paper Authors 

0G-3A 172 Maehl, Israel, Klarmann (WU) 

043 Webber, Lezniak, Kish (UNH) 
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The isotopic composition of Fe and other heavy elements can provide impor- 
tant information on the nucleosynthesis associated with cosmic ray source 
regions. The relative abundances of 54 Fe and 56 Fe are extremely sensitive 
to such details of the source region as temperature, density, and neutron 
excess (see, e.g., Reames, 1970, and Arnett et al., 1971). In addition, the 
production of s °Fe may require source regions different from those producing 
^ 4 Fe and 56 Fe (Truran, 1972). Thus, the isotopic composition of Fe and other 
heavy elements is of great interest. 

The two measurements of Fe isotopes reported at this conference involve quite 
different techniques. The WU group used the geomagnetic cutoff technique 
described by Lund, Rasmussen, and Peters (1971) at the Hobart Conference. 

This method is based on the comparison of the ratio of the flux above a 
geomagnetic cutoff to the flux above a given velocity for one element to the 
ratio for a second element. Since the ratio for a given element depends 
on the spectrum and on the average A/Z for that element, the ratios for two 
different elements will differ only to the extent that either the spectra 
or the average A/Z for the two are different. Using this technique, the 
WU group reports that the average atomic weight for Fe is 

A = 54.6 ±0.5 

The UNH group used an entirely new 
technique which is based on simul- 
taneous measurements of the Cerenkov 
radiation (C) and the total energy 
(E) for each incident particle. This 
technique yields high mass resolution 
over a restricted energy interval 
just above the threshold of the 
Cerenkov radiator. An example of 
the observed mass distributions is 
shown in Figure 3 for a number of 
elements including Al and Fe. Note 
that the mass distribution at Fe is 
rather broad. The authors have 
grouped the observed events as 
follows: 

S3+54 Fe . 15 

55+56+57p e . 23 

58+59+60p e ; 15 


There is no doubt that this Fe distribution and the WU measurement of the mean 
atomic weight are inconsistent. In discussions during and after the sessions, 
this discrepancy was not resolved, although there were comments about the 


Fig. 3. The observed mass 

distributions reported 
by the New Hampshire 
group . 
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statistical weight of the observations and about the background in the 
relatively thick UNH total-energy detector due to nuclear interactions.* 
The interest in these results certainly warrants additional analysis and 
a continuing improvement of these and other techniques. 

I should mention that both of these groups also reported results for other 
nuclei with 10lZ<26. 

V. Cross Sections and Model Calculations. The relevant papers can be 
grouped as follows : 


Session 

Paper 

Authors 

OG-10 

444 

Jung, Suren, Sakamoto, Jacquot, Kaiser 
Schmitt 

OG-11 

406 

Biswas, Ramadurai 

0G-3A 

225 

Tsao, Shapiro, Silberberg 

0G-3A 

- 

Damle, Webber, Kish 

OG-10 

252 

Raisbeck, Perron, Toussaint, Yiou 

OG-10 

253 

Raisbeck, Yiou 


251 

Yiou, Raisbeck, Perron, Fontee 

OG-10 

325 

Silberberg, Tsao 


198 

Ayres i Schmitt, Merker, Shen 

OG-10 

482 

Casse 


The papers in the first group contain calculations of the expected inter- 
stellar abundances of D and ^He. Jung et al. use their recent studies of 
p,a inelastic processes to calculate the interstellar spectra, which they 
modulate using the force-field approximation so that comparison can be made 
with data obtained at 1 AU. Biswas and Ramadurai do a similar calculation, 
but assume a unique source spectrum which they refer to as the Fermi spectrum. 
These authors use a numerical solution of the full modulation equation in 
order to compare with data at 1 AU. It will be interesting to compare in 
detail these and other calculations with each other and with the new D 
results discussed above. 

The second group of papers contain calculations of the secondary production 
of the heavier isotopes. Tsao et al. have continued their calculations of 
abundances and now present results for essentially all of the isotopes from 
Li through Ni. 


*For later comments on this discrepancy, see the Post-Conference Comments 
at the end of this paper . 
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Figure 4 illustrates some of their 
results compared to universal 
abundances . Such results depend on 
many assumptions , many of which can 
be tested by the comparison of such 
detailed calculated abundances with 
isotope data. 

The paper by Damle et al. contains 
calculations of a more restricted 
set of isotopes, specifically those 
of Al, K, Cr, and Mn. These iso- 
topes include 26 Al, 53 Mn, and 5t+ Mn 
which may be useful cosmic ray 
clocks. The two calculations are 
in qualitative agreement on which 
isotopes should dominate each 
element. However, the actual 
fractional abundances of isotopes 
of a given element are different 
in the two calculations, with an 
rms difference of ~ 0.06. At 
present, cosmic ray measurements 
are subject to larger uncertain- 
ities, but this situation could 
change dramatically in the next 
two years. 

The paper by Raisbeck et al, concentrates on the importance of electron 
capture isotopes in cosmic ray studies. The authors discuss the general 
aspects of the problem and the need for specific data such as electron 
attachment and loss cross sections at high energies. They also present a 
calculation of the fraction of surviving as a function of energy. A 
substantial energy dependence is predicted in the 100 to 500 MeV/nucleon 
interval. 

The third group of papers by the Orsay group contains some new results and 
a summary of their current program for measuring the relevant nuclear cross 
sections. Such experimental results will become increasingly important in 
the next few years as the quality and quantity of cosmic ray isotope data 
improve. 

In the absence of a complete set of measured cross sections , we will have 
to rely upon calculations. The fourth group of papers includes results 
of two types of Monte Carlo calculations which may provide more accurate 
predictions for some reaction cross sections. At present, it seems that 
the Monte Carlo calculations are not sufficiently more accurate than the 
semi -empirical calculations to warrant the much greater computation time. 
However , continued development of the Monte, Carlo technique may results in 
significant improvements. 



Moss Number 


Fig, 4 . An example of the 

calculated comic ray 
i8otope abundances 
compared with universal 
abundances (Tsao, et al.). 
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In the last group I have included Casse 1 s calculation of the half life 
of against g + decay. Normally, 56 Ni decay is dominated by electron 

capture, ‘so there is no laboratory data. Casse^ calculates that the g* 
decay half life is - 2 x 10^ y. Thus, if some of the ^Ni thought to be 
present in the outer layers of a supernovae can be directly accelerated 
to high energies, it might be observable in cosmic rays at energies high 
enough so that the time-dilated half life is comparable with the age of 
the source. 

VI. Instrumental Considerations. In addition to some of the above 
observational papers , the following papers from the techniques sessions 
are relevant: 


Session 

Paper 

Author 


TI-4 

435 

Verma, Herzo 

(LSU) 


042 

Webber, Lezniak, Kish 

(UNH) 


346 

Fisher, Ormes, Hagen 

(GSFC-2) 


481 

Casse et al. 



490 

Linney et al. \ 

(SACLAY/DSR1) 


714 

Valot et al. J 



Verma and Her zo describe a new instrument based on the dE/dx^Range-Energy 
technique using a spark chamber for determining the particle path. The 
instrument was flown in summer of 1973 and is expected to have a mass 
resolution of a ~ 0.5 amu for the isotopes of H through Be. 

VT 

Webber et al. provide further details of their hew Cerenkov - total energy 
(C-E) technique which I have already discussed. Presently the charge 
resolution varies from a - 0.4 to ~ 0.8 amu, depending on Z. The authors 
feel a substantial improvement is possible. 

Fischer et al. described a new instrument combining both multiple dE/dx- 
Range-Energy and C-E techniques. The instrument was flown in 1973 and is 
expected to yield a mass resolution of o ~ 0.4 amu for low-Z nuclei and 
~ 0.7 amu for Fe. 

V" v V V 

Casse et al. present details of a new Cerenkov-Rang e-Cerenkov (C-R-C) 
technique which has been calibrated at the Bevatron. Mass resolution from 
the calibrations was a ~ 0.35 to 0.6 amu for N and C. The authors expect 
that the resolution can be further improved. 

Linney et al. describe the new Cerenkov radiators with refractive indices 
in the range from 1.05£n£1.2. The radiators, which are made of compressed 
silica powder, were calibrated at the Bevatron. Such low index radiators 
are a key design element in future isotope studies at high energies. 
Additional calibration details for powder counters as well as for more 
conventional glass and liquid Cerenkov radiators are discussed in Valot et 
al . 



MASS RESOLUTION (cr) (omul 
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In. order to summarize the current status of the measuring techniques for 
isotope studies, I thought it would be interesting to plot the mass 
resolution achieved or predicted for various techniques as a function of 
Z, the charge of the cosmic ray nucleus. Figure 5 contains such a 
summary for "present" and "future" results. The "present" resolution 



Fig. 5, A summary of the mass resolution (a) re-ported for various elements 
(Z) compared with the resolution required to separate adjacent 
isotopes (bn = V with various relative abundances. The " present " 
resolution values are based on cosmic ray data t while the "future" 
resolution values are based on calibrations and calculations . 


results are based on measurements of cosmic ray isotopes as reported in 
the various papers which I have discussed in previous sections. The 
"future" resolution results are based on calculations or calibrations 
described In the instrumentation papers discussed earlier and not on actual 
flight data. I have not attempted to include other "future" techniques 
which were not reported at this conference. 




Figure 5 also Includes 5 resolution levels with labels of 1:1, 2:1, 10:1, 
100:1, and 1000:1. These levels are indicative of the resolution required 
to separate adjacent isotopes (Am=l) with relative abundances of 1:1, 2:1, 
etc. For isotopes separated by 2 amu, the levels would be a factor of 2 
x higher. Thus, separation of two even isotopes with relative abundances 
of 10:1 requires a resolution of a ~ 0.6 amu, provided that the inter- 
mediate odd isotope is non-existent. 


Examination of Figure 5 indicates that satellite instruments are now 
capable of separating isotopes of only the lightest elements (at relatively 
low energies). The major improvements in the immediate future seem to be 
in balloon instruments for heavier nuclei at higher energies where resolution 
of individual isotopes may become a reality. So, the next two years should 
be exciting. Of course, there has been progress in the last two years. 

At Hobart there was controversy about the elemental abundance of Fe. At 
this conference we have progressed to a controversy about the isotopic 
abundance of Fe. Hopefully, there will be even more stimulating results 
at the next conference. 


POST-CONFERENCE COMMENTS . Since the conference I have given further thought 
to the use of the geomagnetic cutoff to determine the average A/Z of dif- 
ferent elements. The actual formulation used by the WU group (paper 172) 
is based on the assumption that the integral momentum-per-nucleon spectrum 
for element i can be represented by ** p"^i over the interval between 

the cutoff and a fixed momentum/nucleon p. It is not possible, of course, 
to determine the correct power law in this interval with a balloon flight at 
a 3 GV cutoff because of the unknown penumbral and east-west effects. 
Therefore, the WU group used the spectral index in a slightly higher energy 
interval, i.e., the - 830 to - 1270 MeV/nucleon interval, as an estimate 
for the index in the 700 to 860 MeV/nucleon interval near cutoff. It 
should be noted that if the absolute value of the actual Yp e in the ~ 700 

to - 860 MeV/nucleon interval was smaller than in the measured energy interval 
then the derived value of A/Z would be too small. 


A qualitative evaluation of this possibility can be obtained by examining 
the differential kinetic energy spectra shown in Figure 4 of paper 171. 

It appears that the shape of the Fe and the 221Z524 spectra are changing 
rather rapidly in the 800 to 1300 MeV/nucleon interval, making it difficult 
to accurately extrapolate into the crucial 700 to 860 MeV/nucleon interval. 
Since the trend in Figure f is toward a flatter Fe spectrum at lower energies, 
i.e., a smaller absolute value for Y Fe> the mean A/Z for Fe may be larger 

than reported. Because of this uncertainty, the WU group has subsequently 
withdrawn their result for the mean mass of Fe. 
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The cutoff technique relies upon accurate knowledge of the spectra In the 
cutoff region, suggesting that essentially simultaneous flights at low 
and high latitudes will be necessary. In this way, the mean A/Z can be 
determined at the higher cutoff using the accurately measured spectra 
obtained at the lower cutoff to specify the y^ in the critical energy 
interval associated with the higher cutoff region. 
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